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A photometric integrating sphere containing a single disk-shaped baffle and a point source with 
an arbitrary angular intensity distribution is described by an integral equation. Numerical solutions 
of this equation are presented for several different source distributions. The errors involved in com- 
paring the various sources are evaluated as a function of baffle size and position. 
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1. Introduction 

The integrating sphere is used in photometry to 
measure total luminous flux. In a sphere with a dif- 
fusely reflecting coating, the illuminance by reflected 
light is approximately the same at all points on the 
sphere surface and is proportional to the total flux 
emitted by the lamp. This reflected light is measured 
at a point on the sphere surface which is shielded 
from the direct rays of the lamp by a baffle. Its in- 
tensity depends somewhat on the geometry and re- 
flectance of the objects inside the sphere and on the 
angular intensity distribution of the light from the lamp. 
Unless this dependence is known, errors will occur 
in measuring the luminous flux from different types of 
lamps. 

It is possible to describe the nonempty integrating 
sphere by an integral equation [l], 1 whose solution 
gives the illuminance by reflected light at every point 
on all the surfaces inside the sphere. The present 
paper describes a numerical solution of this equation 
for a sphere containing a point source with an arbitrary 
angular intensity distribution and a disk-shaped baffle 
located perpendicular to the polar axis of the sphere. 
The existence of an °°-fold axis of symmetry simplified 
the work considerably by making it possible to do the 
azimuthal integrations analytically. The illuminance of 
the surface inside the sphere was calculated for 
several values of the baffle radius and its position along 
the polar axis. Although hemispherical source distribu- 
tions were used for most of the calculations, the results 
are applicable to any small lamp whose intensity distri- 
bution is known. The extension of the technique to 



1 Figures in brackets indicate the literature references on page. 



larger sources, like fluorescent lamps, should not be 
difficult. For convenience, uniform reflectance was 
assumed for all surfaces. The model, however, re- 
quires only that the reflection coefficient have axial 
symmetry about the polar axis. No consideration was 
given here to the effects of nondiffuse reflection. 

2. Theory 

2.1. General Integral Equation 

Let da be an area element of the baffle or inside 
sphere surface about a point labeled a, and let da' 
be another element of area about a point a . Let x 
be the vector directed from a to a, and let n' and n be 
the respective unit outward normals to da and da. 
Let E(a) represent the illuminance at point a by 
reflected light only, and let E Q (a) be the illuminance 
caused by direct light from the source. This latter 
quantity is determined by the overall geometry and 
angular intensity distribution of the source. All surfaces 
are assumed to be perfect diffuse reflectors. Let 
p(a) be the total reflectance at point a. The flux 
dF(a, a') coming from da' to da is given by the 
expression 

dF(a,a')=p(a') [E (a) +E(a')] 

da (-x-n) da (x • n) X S(a, a')l (ttx 4 ). (1) 

A derivation of this equation is given by Jacquez and 
Kuppenheim [2]. In the present work a factor S(a, a') 
has been introduced to account for the screening of 
one part of the sphere's surface from another by the 
baffle. If all points on da can be seen from a, then 
S(a, a) = 1. If da' is completely or partially screened 
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from a by the baffle, then S (a, a') will be zero or some 
number less than unity. The illuminance £ (a) of da by 
reflected light is then equal to the integral of the 
quantity dF{a, a') Ida over all the area elements 
da' on the inside of the sphere. This, of course, in- 
cludes the baffle surface as well as the sphere surface. 
The unknown quantity E(a) is thus expressed as an 
integral containing itself and quantities which are 
known. We have then, 



E(a)= f p(a')[Eo(a') 

J a' 



+ E(a')]T(a,a')S(a,a')da'l7T (2) 



vhere 



T(a,a') = (-x-n')( X -n)/* 4 



(3) 



The numerical solution of (2) can be obtained by 
first dividing the surfaces involved into a number of 
small area elements 8a' having roughly the same size, 
and approximating the integral by a sum over these 8a'. 
Values of the functions S(a, a') and T(a, a') for all 
possible pairs of area elements are fixed by the geom- 
etry of the system, while p(a) and Eo(a) are assumed 
to be specified. One then guesses a solution for E(a), 
say E°(a), and uses (2) to get a first approximation 
E 1 (a). The process is repeated using this first approxi- 
mation in (2). Repetition of this procedure gives a se- 
quence of functions E l (a) which may converge. If it 
does, the accuracy of the solution may be checked by 
repeating the whole calculation with smaller area ele- 
ments 8a\ or with a better approximation to the inte- 
gral, to see whether the limiting value of the function 
E(a) changes significantly. The former method was 
used in the present work. 

Equation (2) is quite general except for the restriction 
that the surfaces be diffuse reflectors. In practice, the 
difficulties involved in evaluating S (a, a') and T(a, a') 
make it desirable to keep the geometry as simple as 
possible. This will also reduce the number of terms 
needed in the representation of (2) as a sum, thereby 
saving time on the iterative calculations. In particular, 
the use of an experimental arrangement having an 
°°-fold (C x ) axis of symmetry greatly simplifies all as- 
pects of the analysis. Furthermore, this symmetry re- 
striction does not impose any symmetry requirements 
on the angular intensity distribution of the source. It 
need be applied only to the geometry and reflectance 
of the surfaces in the sphere. 

2.2. Specific Model Used in the Calculation 

a. Classification of Surfaces 

The sphere geometry in the present calculation was 
chosen to approximate the substitution method of using 
the sphere. Positions of the various surface elements 
were specified in spherical coordinates with the origin 
at the sphere center. This model had the following 
characteristics: 



(a) A point source at the center of the sphere of ra- 
dius /?, with no restrictions on its angular intensity 
distribution. 

(b) A flat, infinitely thin, circular baffle of radius B 
with its center on the polar axis a distance D from the 
sphere center and oriented normal to the polar axis. 

(c) Uniform reflectance for sphere and baffle 
surfaces. 

For this model, the solution of (2) requires the evalua- 
tion of five integrals. Let the symbol s designate points 
on the sphere, «a points on the north side of the baffle, 
and as points on the south side. Let dE = dE(a, a') be 
the illuminance at point a produced by the area ele- 
ment da' around a'. Equation (2) then takes the form 

E(s)= J dE(s, s') + I dE(s, a N >)+ J dE(s, a s >) 

Js' J a N t J a S' 

E(a N ) = dE(a N , s') 

Js' 



E(a s ) 



-I 



dE(a s , s') 



(4) 



Numerical integration of (4) is necessary only for 
the polar angle 6 since the integration over the azi- 
muthal angle <p can be done analytically. Throughout 
the discussion, points on the sphere surface will be 
specified by the coordinates 6 and <p written without 
subscripts. The radial coordinate r is a constant equal 
to R for these points. Baffle points will be designated by 
the subscript B on the spherical coordinates. The 
radial coordinate r B is related to 6b by the equation 



r B = D cos 9 B . 

b. Separation of the Azimuthal and Polar Integrations 



(5) 



When both a and a' are on the sphere, T(a, a') is a 
constant. If one of the points is on the baffle, the de- 
termination of S(a, a') is trivial. To carry out the inte- 
gration involving the azimuthal variable it is therefore 
only necessary to calculate the azimuthal dependence 
of S(a, a') when both points are on the sphere, and 
that of T(a, a') when one point is on the baffle. These 
calculations are outlined in appendices I and II. Fol- 
lowing this, Eo(a) for the different surfaces must be 
calculated in terms of the angular intensity distribution 
of the source. This is done in appendix III. For con- 
venience, it will be assumed that the source intensity 
distribution is independent of the azimuthal angle <p. 
This means that E(a) will also be independent of cp. It 
is important to note, however, that the value of E(a) at 
a pole of the sphere depends only on the values of Eo(a) 
averaged over <p, provided the reflection coefficient has 
C x symmetry. If the observation point is at a pole there 
are no restrictions on the azimuthal variations in the 
source intensity. This point will be discussed in ap- 
pendix IV. 

c. Method of Integrating Equation (2) over the Polar Angle $ 

Because of the C x symmetry of the sphere and baffle, 
it was possible to perform the integration of (2) over 
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the azimuthal angle <p analytically. This integration is 
discussed in appendices II and III. After this integra- 
tion the five integrals in (4) become 

f dE(s, s')=[ ' p(0')[E o (0') 

Js' JO 

I 

+ E(0')]a(0 9 6') sin 0'd0'l{4m) 
J dE(s, a N >) = \ S pn(0b')E n (0 b >)t(0, 6 b , ) D 2 

Ja N , Jo 

k 

cos 6 B > sin 6 B >d6 B >l7T, 6 < 6 e 

\ 

- I dE(s,a S ')= \ 6s Ps(Ob>)[Eo(Ob>)+E s (O b ,)] 

J a S' Jo 

r(0, B >) X D 2 cos B > sin 6 B >dO B >liT, 6>6 e 
[dE{a N ,s')= V' p{d')[Eo{6') 

j* +E($')] t(0\ 6 b ) R 2 sin 8' d6' I tt 



" I dE(a s ,s')= [ n p(8')[E (0') 
Js' Je p 



+ £(©')] t(0\ 6 b )R 2 sin O'de'l-rr. 



(6) 



The functions E (6), E (6 B ), cr(0, (9'), and r(0, fl ) 
have been denned by (A7), (A8), (A3), and (A5), re- 
spectively. Reflection coefficients are designated by 
p(0), Pn{0), and ps(0) for sphere points, north 
baffle points, and south baffle points. The functions 
2?(0), En(Ob), and E S (6 B ) are averages over the 
azimuthal angle of the indirect illuminance for the 
sphere and baffle points. The parameters 6 S and 6 e are 
discussed below. The next step, integration over the 
polar angle, was done numerically, by a modified form 
of the trapazoidal rule. 

To do this, the range covered by 6 was divided into 
20 (or 40) equal intervals 86. These are shown in figure 1 
along with the other sphere and baffle parameters, 
D, B, R, 6 S , d e , and 6u The angle S has the value 
tan -1 (Z?/Z)); it defines the portion of the sphere 
shielded from the direct rays of the source. Thus, 
we have Eo(0)=O for 6<6 S . The angle 6 e is deter- 
mined by the intersection of the plane of the baffle 
with the surface of the sphere. It has the value cos -1 
(D/R). Sphere points for which 6<6 e cannot be 
seen from the south side of the baffle, while those 
for which 8 > 9 e cannot be seen from the north side. 
We have adopted the restriction that 6 e must equal a 
value of 6 which marks a boundary between two 
86 intervals (see figure 1). This in turn limits the 
number of D values that may be used. This is an un- 
necessary restriction which could be relaxed. The 
angles Si are defined by the equation 




FIGURE 1. Definition of sphere parameters. 

R is the radius of the sphere; B is the baffle radius, and D its distance from the center of the 
sphere where the light source is located. The positive z axi9 passes through the north pole. 
The range of the polar angle 6 is divided into 20 equal intervals 80. Their centers are denoted 
by the polar angles 0i, i = 1 . . . 20. 6, defines the portion of the sphere shielded from the 
direct rays of the source. 6 e is determined by the intersection of the plane of the baffle with 
the surface of the sphere. 

They thus lie at the midpoints of the 86 intervals. In 
figure 1, the circles on the perimeter of the semicircle 
indicate the Si values, and the other index lines show 
the boundaries of the 86 intervals. 

The range covered by 6 B was not divided into equal 
intervals. Instead, the baffle radius was divided into 4 
equal parts, of length 8B = BI4>. The centers of these 
intervals fall at distances Bi from the center of the 
baffle. The Bi are given by the equation Bi = (i — i)8B, 
i = 1 . . . 4. The values of 6 B i corresponding to these 
positions are given by 



Bi = tan- 1 (BJD) 



i = 1 . . . 4. 



The increment 86 B thus depends on i; 



86 b 



8B(D cos 2 6 Bi ) 



1 



4. 



(8) 



(9) 



6i=(i-i)88 



i=l 



20(or40). (7) 



Before representing the integrals in (6) by summa- 
tions, let us examine how rapidly the different functions 
in the integrals vary over intervals of size 86. A study 
of figure 1 will show that the functions E(6), En{6 b ), 
Es(6 B ), and t(0, Ob) experience only small relative 
changes over any of the 86 or 86 B intervals. The func- 
tions E (6), E (6 B ), and cr(6, 6'), however, can 
change markedly over certain 86 intervals. For instance, 
E (6) has a discontinuity in the 86 interval which con- 
tains 0.s-. To cope with these rapid variations, the slowly 
varying functions were assumed to be constant over 
the interval 88 and were taken out from under the inte- 
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gral sign. The remaining portion of the integral was 
then approximated by dividing 86 into 10 equal parts 
and summing over these smaller intervals. The same 
procedure could be used for the intervals 86 B if ever 
there is a large variation in E ( 6 B ) . 

Written as a summation, the integral equation (4) has 
the form 



m 



Ei = ^[(DAA) ij +(AA) ij E j ] 
j=i 

+ 2l(AN)vENj+(DAS) v 

+ (AS)ijESj] i=l . . . 
EN l = f i [(DNA) iJ + (NA)tjEj] 

m 

ES i =^[(DSA) ij +(SA) ij E j ] 

j=e+l 



m 



i= 1 



i= 1 



. 4 



(10) 

where Ei is the value of E(6) when 6= Si, ENi and 
ESi are the values of En (6b) and Es (6b) when0# = 0#j, 
e = 6 e j8d, m = 20 or 40, and the various coefficients 
are defined as follows: 

(AA) ij = p(6 j ) j o-(0 u 6) sm6d6/(47r) (11) 

where j±=6j±i86, and the integral is evaluated 
numerically; 

(AN)ij = p N (6 Bj )r(6 u 6 Bj )D 2 cos 6 Bj sin 6 Bj 8d Bj l7r 

i = 1 . . . e 

= i = e+l . . . m (12) 

(AS)ij = Ps(6bj)t(0 u 6 Bj )D* cos 6 Bj sin Bj 80 Bj lir 

i = e+1 . . . m 

= i = 1 . . . e (13) 

(NA)ij = p(6j) T (6 h 6 Bi )R 2 sin 0,-80/tt 

J=l . . . e 

= ;= e+1 ... to (14) 

(S^)y = p(6j)T(6 h 6 Bi )R 2 sin 0,-80/tt 

jf = e + 1 . . . m 

= j= 1 ... e. (15) 

The coefficient (DAA)ij is given by (11) with the addi- 
tion of Eo(6) under the integral sign. We can write 
(DAS)ij = E (6 B j)(AS)ij unless E (6 B ) shows appreci- 



able variation over 50#j, in which case (DAS)ij should 
be written as an integral over 86 B j. The remaining 
coefficients (DNA)ij and (DSA)ij are equivalent to 
the expressions (14) and (15), each multiplied by£' (^) 
and integrated over the ran gey' ± as in (11). 



2.3. Solution of the Integral Equation 

As outlined in (IIA), an iterative method was used 
to solve (10), which is the integral equation (2) ap- 
proximated by a summation. Equation (10) contains 
the unknown quantities Ei, ENi, and ESi and a series 
of coefficients which depend on the geometry and 
reflectance of the surfaces, and on the angular 
intensity distribution of the source. 

A convenient zeroth order approximation is the 
solution of (2) for a sphere containing only a point 
source with a uniform intensity distribution. For area 
elements on a spherical surface, T(a, a') = l/(4/? 2 ) 
for all a, a'. Without the baffle, we also have S(a, a) 
= 1 for all a, a . The integral (2) in this case becomes 

E(a) = \ p(a') [Eo(a') +E(a')} da'K^R 2 ) 

Ja' 

(16) 

where the integration is over the sphere surface. Since 
this integral contains nothing which depends on a, 
we must have E(a) independent of position. Letting 
E(a) = £, we can solve (16) for E to get 

E= [ p(a')Eo(a')da'\l- [ p(a') da'/^TrR 2 )] 

(47T/? 2 )- 1 . (17) 

For uniform reflectance p, and a centrally placed 
source with a uniform angular intensity distribution, 
one has Eo(a) =Eo, and (17) becomes 



E = pE /(l-p). 



(18) 



The value E calculated from this equation was used 
as the first approximation to all the unknowns, £"i, 
EN u and ESi. 

With p = 0.9, values of E(a) satisfied (10) to better 
than 1 part in 10 5 after about 35 iterations on a high 
speed computer for both the 24 and 44 point models. 
The solutions were checked to see that they satisfied 
the condition of flux balance. This is given by the 
relation 

f E (a) da= [ [1 -p(a)][Eo(a) +E(a)]da 



where the integrations are over all the surfaces in the 
sphere. The integral on the left is the total flux into 
the sphere, and that on the right, the total flux out. 
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3. Calculations and Results 

3.1. Description of Source Distributions 

Two source intensity distributions, each uniform in 
one hemisphere and zero in the other, were used for 
most of the calculations. Specifically, these distribu- 
tions were 

I N (6) = 2 lm/sr ^ 6 < tt/2 

= IT/2 ^ 6 ^ 77 

I s (6) = ^ 6 < tt/2 

= 2 tt/2 ^ 6 ^ 7T. (19) 

The total flux F<> in each case is, by (A10), equal to 
47rlm. 

With the sphere geometry discussed in the next 
section, an additional group of distributions having 
the form 
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FIGURE 2. Variation with the polar angle of the indirect illumi- 
nance K{6)for the northward directed hemispherical source distri- 
bution, and the reference sphere geometry. 

E {d) is the direct illuminance. 



= 



(all other 6 values) (20) 



was used. The ft are given by (7). The total flux is 
the same for each distribution and is equal to 47T lm. 
We shall discuss later how the results obtained from 
this series of distributions can be used to determine 
E(a) for an arbitrary source distribution. 

3.2. Results for a Reference Sphere 

For much of the work, a particular set of sphere 
parameters was used. They were chosen to correspond 
roughly to those of a typical experimental sphere. 
They were, R = 1, B = 0.1, and D — 0.31 m. A uniform 
reflection coefficient of 0.9 was used for all surfaces. 
Calculated values of the indirect illuminance E(6) 
for the 24-point model are shown as polar plots in 
figures 2 and 3 for the north and south hemispherical 
source distributions, respectively. If there were no 
baffle, E(6) would have a uniform value of 9 lm/m 2 . 
The E{6) values shown in figures 2 and 3 are lower 
than this on the average because of light absorption 
by the baffle. The E(6) values are enhanced on the 
side of the sphere opposite that being directly illumi- 
nated. Values of the indirect illuminance of the baffle 
show the same effect. They are shown in table 1 for 
the reference sphere and the hemispherical distribu- 
tions. For the north distribution, the north side of the 
baffle, being opposite the directly illuminated surfaces, 
receives more indirect illumination than the south side. 
The opposite situation holds for the south distribution. 
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FIGURE 3. Variation with the polar angle of the indirect illumi- 
nance K(6)for the southward directed hemispherical source distri- 
bution, and the reference sphere geometry. 

E (0) is the direct illuminance. 
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Table 1. Indirect illuminance of baffle for reference sphere 
dimensions, p = 0.9, hemispherical source distributions, and 
24-point model 



/.v(0)(lm/sr) 


/ 5 (0)(lm/sr) 


i 


<Mdeg) 


ENi(\mlm 2 ) 


ESi(\m/m 2 ) 


ENi{\m/m 2 ) 


ESi(lmlm 2 ) 


1 
2 
3 
4 


2.32 

6.94 

11.46 

15.86 


9.5510 
9.5522 
9.5547 
9.5583 


8.2733 
8.2736 
8.2743 
8.2754 


8.0889 
8.0888 
8.0886 
8.0883 


9.7416 
9.7413 
9.7407 
9.7398 



This table also shows the large difference in indirect 
illuminance of the north and south sides of the baffle, 
and its large dependence on the source distribution. 
Actually, it is the large change in ENi which is pri- 
marily responsible for the change in E\ ~ E(6 = 0) 
(see figs. 2 and 3) when the distribution shifts from one 
hemisphere to the other. The changes in the amount of 
light coming to the 6\ position from the other surfaces 
are small by comparison. 

Because the detector would normally be at the north 
pole, it is the change in Ei that is of primary interest. 
For the reference sphere, E\ changes by 0.25 percent 
when the hemispherical source distribution is switched 
from north to south; values of E(6) at other 6 values 
change to a greater or lesser degree. 

3.3. Effect of Changing the Baffle Position 

The effect of changing the position of the baffle along 
the polar axis is shown in figure 4. Here, values of 
E(6) are plotted as a function of 6 for four different 
baffle positions D, including that occurring in the 
reference sphere. The radius of the baffle was held 
constant at Z? = 0.1 m. For each case, there are two 
curves corresponding to the two hemispherical dis- 
tributions. That arising from the north distribution is 
labeled N. Of the four positions, the difference in E\ 
between the two source distributions is smallest for 
the baffle position used in the reference sphere. This 
position does not, however, minimize this difference. 
A plot of Ei(N) — Ei(S) versus D at constant B, where 
N and S designate the distributions, is shown in 
figure 5 by the filled-in circles. This plot shows that 
for D ~ 0.28 m there would be no difference at all in 
Ei between the two distributions. This optimum baffle 
position differs considerably from the value of D = 
0.372 R, which minimizes the screened areas of the 
sphere surface [3]. Minimizing the screened areas does, 
however, minimize the change in E{6) averaged over 
the sphere that occurs when the hemispherical source 
is directed from north to south. The sum of the squares 
of the differences Ei{N)—Ei(S), each weighted by 
sin 0i to account for the different sizes of the sphere 
areas represented by the 0j, was smaller for Z>=0.38 
than for D = 0.31, 0.45, or 0.28 m. (Values ofE0) and 
Ei(S) for D=0.28 and 0.38 m were estimated by 
interpolation between the curves shown in fig. 4.) Thus, 
while the dependence of the E(6) values around the 
sphere on the source distribution is minimized by 
D = 0.372, the experimentally important quantity 
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Figure 4. Indirect illuminance E(0) versus 6 arising from the north 
{labeled N) and south hemispherical distributions for different 
values of the baffle position. 

The other parameters have reference sphere values. 



£i(/V)-£i(S) is minimized by D=0.2S m. At this 
stage in the calculations it is not possible to generalize 
this result to other source distributions. The optimum 
baffle position may not be the same for source distribu- 
tions differing from the hemispherical ones considered 
here. 



3.4. Effect of Changing the Baffle Radius 

At a fixed value of Z), decreasing the baffle radius 
will decrease the error Ei(N) — £\(S) made in 
observing the hemispherical distributions. Examples 
of this are shown in figure 5. The three points con- 
nected by the dotted line correspond to combinations 
of B and D which keep constant the area of the sphere 
screened by the baffle from the observation port. 
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FIGURE 5. Difference Ei(N)— Ei(S) between the indirect illumi- 
nances near the north pole produced by the north and south hemi- 
spherical distributions as a function of baffle radius and position. 



3.5. Determination of E(H 0) for an Arbitrary 
Source Intensity Distribution 1(6) 

Calculations of E\ using the reference geometry 
with the 24-point model were done for a number of 
the distributions U(6) defined by (20). Several of the 
E\ values were obtained from the others by graphical 
interpolation. All are shown in table 2. 

Consider a source with an intensity distribution 
1(6) which can be approximated by a linear combina- 
tion of the Ii(6) distributions. Let 

20 

/(«)=y»i/i((i), 

where the weights w i satisfy the relation 

20 
i=\ 

Then, if we define E\(i) as the value of E x calculated 
for the distribution Ii(0) we have simply 



20 



WiE x (i), 



where E\ is the value for the distribution 1(6). Values 
of E\(i) from table 2 were used with appropriate 
weights to calculate E\ for the hemispherical dis- 
tributions In(6) and Is{6). The results differed by 
0.03 percent from those calculated directly from In{6) 



and Is{6). The use of interpolated values for some of 
the E\(i) is probably responsible for most of the 
difference. 



Table 2. Indirect illuminance of north pole region for reference 
sphere dimensions, p = 0.9, source distributions Ij(#), and 24-point 
model 



i 


Bi 


IiW 


Ei 




(deg) 


(lm/sr) 


(lm/m 2 ) 


1 


4.5 


159.4 


7.8329 


2 


13.5 


54.65 


8.1143 


3 


22.5 


33.33 


9.0464 


4 


31.5 


24.41 


9.0216 


5 


40.5 


19.64 


9.0010 


6 


49.5 


16.77 


8.9823 


7 


58.5 


14.96 


8.9666 


8 


67.5 


13.81 


8.9542 


9 


76.5 


13.12 


a 8.951 


10 


85.5 


12.79 


a 8.950 


11 


94.5 


12.79 


8.9500 


12 


103.5 


13.12 


a 8.950 


13 


112.5 


13.81 


a 8.949 


14 


121.5 


14.96 


8.9487 


15 


130.5 


16.77 


a 8.947 


16 


139.5 


19.64 


a 8.944 


17 


148.5 


24.41 


8.9387 


18 


157.5 


33.33 


8.8056 


19 


166.5 


54.65 


8.4642 


20 


175.5 


159.4 


8.0655 



a These values were obtained by graphical inter- 
polation between the unmarked values. 



4. Discussion 

The integral equation description of the integrating 
sphere offers an accurate way to determine the errors 
caused by the screening effects of objects inside the 
sphere. While its numerical solution is possible, in 
principle, for any sphere geometry, a great simpli- 
fication results when an °°-fold axis of symmetry exists 
for the various surfaces and their reflection coefficients. 
The computation time is so short for such a system, 
that it becomes feasible to improve the accuracy by 
using smaller area elements and to extend the calcu- 
lations to spheres containing several objects. 

The effect of increasing the number of sphere points 
from 20 to 40 is shown in table 3 for the reference 
sphere dimensions. The biggest change in the Ei 
values occurred forfi. There were decreases in E\(N) 
and Ei(S) of around 0.03 percent and an increase in 
their difference E\(N) —Ei(S) of 3.1 percent. Further 
improvements in the accuracy of the calculation could 
be achieved by using more points, and/or by using a 
better approximation to the integral. 

An extension of the method to spherical and cylindri- 
cal sources located anywhere along the sphere axis 
appears feasible. This would greatly increase the 
number of systems to which the method could be 
applied. 
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Table 3. Effect on indirect illuminance near* the north pole of The azimuthal dependence ofS(6; 0', <p') is thus 
decreasing the interval size 86, for reference sphere dimensions, 
p = 0.9, and hemispherical source distributions 



EitfV) 
Ei(S) 

£i(AO-Ei(S) 



20 sphere points 



8.9277 lm/m 2 

8.9087 

0.0190 



40 sphere points b 



8.9259 lm/m 2 

8.9063 

0.0196 



Change 
(percent) 



0.020 
0.027 
-3.1 



a The E(0) values given here are those at = 4.5°. 

b The values at = 4.5° for the 40 sphere point case were obtained 
by linear interpolation between the E{6) values calculated at 
(9 = 2.25 and 6.75°. 



S(0; 0', if') = 1 for ^ <p' ^ if' s and 

2tt-<p' 8 
= for <p' 8 < <p' < 2n - (p' 8 
= for ^ < Be < 0' =i 77. 

For the integral over <p f we have 



<p' ^ 2tt 



cr(0, ^') 



r 

Jo 



S(0;6',<p')d<p' = 2<p' s . (A3) 



The writer thanks S. R. Kraft for suggesting the The azimuthal integration in the evaluation of the 

integral equation approach to the baffle problem and first integral^ f S 'dE(s, s f ) in (6) is thus accomplished 

for guidance in the initial phase of the work. Ac- since S (a, a') is the only factor in the integrand which 

knowledgment is also given to W. B. Fussel for many contains <p . 
valuable discussions. 
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6. Appendix I: Azimuthal dependence of 
S(a, a') when both a and a' are on the 
Sphere Surface 

Consider two points (0, <p) and (0' , <p') on the sur- 
face of the sphere. These points determine a line which 
intersects the plane containing the baffle at some point. 
Let )8 be the distance in the plane of the baffle between 
this point of intersection and the polar axis. Because of 
the Coo symmetry of the system, the azimuthal depend- 
ence of /3 can be specified in terms of the difference 
<p — <p'. One can thus let <p=0 and calculate /3 as a 
function of the three parameters 6,6', and <p' . It can 
be shown that 



7. Appendix II: Azimuthal Dependence of 
T(a, a') When One Point Is on the Baffle 
and the Other Is on the Sphere 

Consider a point (6, <p) on the sphere surface and a 
point (6b, <Pb) on the baffle. It can be shown that 



T(6,<p; d B ,<PB)=k[fi + gi cos (<p — <p B )] 

[/ 2 + #> cos (<P-<Pb)]~ 2 
where 



(A4) 



/3 2 =f+g cos <p' 



(Al) 



where 

/= iR 2 [ (1 - O 2 sin 2 + (1 + T) 2 sin 2 0' ] 

g = i« 2 (l - r) (1 + O sin 6 sin 6' 

r=[2D-R(cos 0+cos 0')P(cos 0'-cos 0)]" 1 . 

The line between (0, <p) and (6\ <p') will intersect the 
baffle under the conditions, 0^ 6 < 6 e < 6' =§= 7T, and 
)3 2 = B 2 , where 6 e is the polar angle determined by the 
intersection of the baffle plane with the surface of the 
sphere (see fig. 1). Let <p' s be that value of <p' for which 
(3 2 = B 2 . From (Al) we have 



*>,' = cos- 1 [(B 2 -f)lg}. 



(A2) 



A; = ±(Z>cos 2 fi -/?cos0) 

fi = R-D cos 2 6 B cos 

g\ = — D cos 6b sin 6b sin 

f 2 = R 2 + D cos 2 6 B (D - 2R cos 0) 

gi = — 2RD cos 6 B sin 6 B sin 0. 

In the expression for A:, the negative sign is used for 
baffle points on the north side, and the positive sign 
for points on the south side. 

As one would expect from the C x symmetry, 7\0, <p; 
6b, <Pb) depends on the difference <p — <Pb between the 
azimuthal angles. The integral over either azimuthal 
angle is 

T27T 

t(6,0 b )= T(6,<p;6 B ,<PB)d<p 
Jo 

C27T 

= T(6,<p;6 B ,<pB)d<p B 
Jo 

= 2Trk(f 1 f 2 -g 1 g 2 ) (fl-gl)- 312 . (A5) 
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This function t(6, 6b) furnishes the azimuthal portion 
of all the integrals in (4) except the first. 

8. Appendix III: Determination of the Direct 
Illuminance E (o) 

Let 1(6, <p) be the angular intensity distribution of 
the point source in lumen/steradian. The direct 
illuminance Eo(6, <p, r) at a point (6, <p, r) is given 
by the expression 



Et>(6,<p, r)=/(0,^)(r-n)r- 3 



(A6) 



where r is the vector from the source (located at the 
origin) to the point and n is the unit outward normal 
from the area element about the point. For sphere 
points (r • n)r -3 has the value l//? 2 , while for baffle 
points it has the value 1/ (D 2 cos ^).When p(a) has 
Coo symmetry, E(a) at the poles is independent of the 
azimuthal dependence of 1(6, (p). If the observation 
point is at a pole it is permissible (see appendix IV) 
to use the average of Eo(6, (p) over <p in the 
calculation. 

For points on the sphere surface, excepting those 
screened from the source by the baffle, we have 



! *> I 



Eo(0) = (1/277-) Eo(0,<p)d<p 



= (U2ir)f*"l(6,<p)dplR*. 



(A7) 
For baffle points, we have 

Eo(Ob) - (l/2ir) J" E (6 B , <p„)d<p B 

= (D 2 cos 6 B ) "' (l/2ir) f ** I(6 B , Vb) d*p B . 

(A8) 

The north, south designation is not needed here, since 
Eo(0 B ) is zero on the north side of the baffle. One 
need only know 1(6, (p) averaged over <p\ i.e., 

1(6) - (II2tt) J*" I(6,<p)d<p. (A9) 

The total flux F from this source is given by the 
integral 



where d£l is the element of solid angle and the integra- 
tion is over \ir steradians. 

9. Appendix IV: On the Absence of Symmetry 
Restrictions on the Source Intensity 
Distribution 

When the reflection coefficient has C x symmetry, 
and the observation point is at the pole, there are no 
symmetry restrictions on the angular intensity dis- 
tribution of the source. To see why this is so consider 
the integral equation (2) for points on the sphere near 
the north pole. 

J'TT [" T27T 

pie 1 ) [E (6',<p') 
o LJo 

+ E(6\ <p')]S(6, ip; 6', <p')dtp f X sin 6'd6'l^7T 

+ f ' Pn(6 b >) \ l J E n (0 B ; <Pb>)T(6, if; 6 f , <p')d<p B ] 

XD 2 cos B > sin 6 B >d6 B >l'7r 
The average of E(6, tp) over <p is 

E(6) = (1/2tt) I 2 " E(6, <p)d<p 

Jo 

= ( p(6 f )[E (6 f )+E(6 f )]a(6,6') sm6'd6'l47r 



f 

Jo 
■ f" S Pn(6 b .)E n (6 B ') 

JO 



Fo= f 1(6, <p)dQ.=-2ir \ 1(6) sin 

J4n J» 



(A10) 



t(0, b -)D 2 cos 6 B r sin B 'dO'lir 

(All) 

where we have used the relations (A3) and (A5). The 
function Eq(6) is the average of E (6, (p) over <p as 
defined by (A7), and En(6b) is defined for the baffle 
points in the same way that E(6) was defined for the 
sphere points. This expression shows that E(d) values 
averaged over the azimuthal angle will satisfy the 
integral equation. 

The flux striking a small round observation port of 
radius p, centered on the north pole of the sphere, is 
given by the integral 

[ 21 \[ dl> E(6,<p) sin 6dd\ dip = 2tt j $IJ E (6) sin 6d6 

where the upper limit of integration is equal to sin" 1 
(p/R). As (All) shows, this integral depends only on 
the average azimuthal dependence of the source 
intensity distribution. 
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